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Abstract 

In this paper, we construct a theory of the NJL-type where superconductivity is 
^ \ present, and yet the super-conducting state remains, in the average, color symmet- 

^ , ■ ric. This shows that the present approach to color superconductivity is consistent 

with color singlet-ness. Indeed, quarks are free in the deconfined phase, but the 
^ , deconfined phase itself is believed to be a color singlet. The usual description of the 

' color superconducting state violates color singlet-ness. On the other hand, the color 

On ■ superconducting state here proposed, is color symmetric in the sense that an arbi- 

trary color rotation leads to an equivalent state, with precisely the same physical 
properties. 



<3\ 
O 

O ■ 1 Introduction 



It is presently accepted that quark and gluon fields are the building blocks of 
hadronic matter, in the framework of quantum chromo dynamics (QCD). The inves- 
ts \ tigation of the phase structure of hadronic matter is a topic of great current interest. 

A diversity of phases is expected at high densities: chiral-symmetry restoration, de- 
confinement and color-superconductivity. Since, due to the complexity of the theory, 
it is extremely difficult, if not impossible, to obtain exact results directly from QCD 
when perturbation theory cannot be applied, effective models, such as the Nambu- 
Jona-Lasinio (NJL) model [U [21 [3], must be employed and have been used with 
great success to investigate the properties of hadronic matter and to develop insight 
into its phase diagram [3]. 

Recently, the color superconducting phase in quark matter has been investigated 
by many authors [5], in the framework of the Bardeen-Cooper-Schrieffer (BCS) 
approach, which is familiar from condensed matter physics. For a recent review, 
see [6]. Quarks are free in the deconfined phase, but the deconfined phase itself is 
expected to be a color singlet. In Refs. [71 [8] it has been argued that in QCD the 
superconducting phase is automatically color symmetrical. Our basic assumption 
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is, therefore, the existence of globahy color-symmetric superconducting phases, and 
our aim is to discuss how these phases may be described in terms of effective models 
with 4 fermion interactions, such as the NJL-like model considered in Refs. [5l [6], 
for which it is assumed that the gluon degrees of freedom have been integrated over. 

A BCS state |$) describes a physical state with zero net color charge if A^i = 
A''2 = N3, where Ni denotes the average number of quarks of color i. This means 
that 

(cJ.|5a3|'J>) = («1>|5a8I^) = 0. (1) 

Here, Sx^ denotes the SU{3) generator associated with the Gell Mann matrix A^. 
However, the requirement ([T]), which is implemented in OITO], is not sufficient to 
insure that |$) is physically acceptable. A stronger condition must then be imposed. 
Indeed, the SU{3) symmetry, being a gauge symmetry, cannot be broken, according 
to the discussion in [71 [8], so that color rotated BCS states must be equivalent in 
the sense of the physics they describe. Let Uc denote an arbitrary color rotation, 
i.e., Uc = exp^^^j^ ix/c^Aj., the parameters Xk being arbitrary and real. The BCS 
state 1$) must be equivalent to the state Uc\^), for any Uc, as far as expectation 
values of physical observables are concerned. Therefore, the condition ([T]) must be 
replaced by 

{m^,sx,Ucm = {^^cSxsUcm = 0, 

for an arbitrary Uc, and this implies 

($|SaJ^) = 0, for A; = 1,2,--- ,8. (2) 

This is the condition the BCS state must satisfy in order to be physically 
meaningful. If only the condition ([1]) is implemented, and not the condition ([2]), the 
BCS state |^>) is, in general, not equivalent to the state Uc\^), so that it describes 
a state belonging to a representation of SU{3) other than the singlet one, which is 
physically unacceptable. In [11] it is shown how the condition ([2]) may be easily 
implemented. 

In the present paper, we apply to the NJL model the new BCS approach devel- 
oped in Ref. [TT] which uses the generalized Bogoliubov transformation and leads to 
a color symmetric BCS vacuum. In Section 2, a mean-field constrained Hamiltonian 
appropriate for the description of color superconductivity is presented. In section 3 
we compare the new superconducting state which satisfies ([2]), with the usual one, 
which is not required to satisfy ([2]). In Section 4, we draw some conclusions. 

2 Color symmetrical superconductivity 

We wish to focus on the diquark condensate {i^^ i^^\jT2il^) , j G {2,5,7}, and, 
for simplicity, we will neglect the quark-antiquark condensate {i/jil^). The notation 
is the usual one, the superscript C denoting charge-conjugation. In this sense, the 
model is not fully realistic, but is adequate for the present development. We assume 
that the pairing interaction is antisymmetric in the color indices and in the iso-spin 
indices, i.e. it involves the Gell-Mann matrices \j, j E {2,5,7} and the iso-spin 
matrix T2 associated with fiavor, when 2 is the number of fiavors. Having in mind 
a Hamiltonian of the NJL [2] type, the mean-field constrained Hamiltonian reads: 
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Kmfa = / d X 



'0(p-7 + M-/i7o)^ + ^ ^ {A*i!^i-f5T2Xj'^ + h.c.) 

j6{2,5,7} 



I A -12 

ie{2,5,7} 



(3) 



where 



A, = -2Gc(V''^75r2A,V'), (4) 

denotes the BCS gap. We use the symbol K, instead of the more usual symbol H 
to stress that this Hamiltonian is constrained in the sense that it fixes the Fermion 
number through the chemical potential fx, which behaves as a Lagrange multiplier. 
Thus the expectation value of Kmfa is the thermodynamical potential which is 
equal to —PV, where P is the pressure and V the volume, and determines the 
equation of state. By {X) we denote the average of X in the BCS vacuum which 
will be specified in the following. The notation is essentially the same as in |12j . 
slightly modified, as is required in order to treat the 3 colors on the same footing. 
In momentum space we have, 

Kmfa = [{Vp"^ + M'^ - ^Cp,,,j> + {Vp'^ + M'^ + A^)Cp,,,j>Cp,r,,i,r 

I P,Vj,k,T,T' 

where Cp ^ ^ ^ and cj, ^ ^ ^ create, respectively, a quark and an antiquark of momen- 
tum p, helicity r], color index j and isospin index r, Cp,?; = — C-p,/y = —Cp-r] = 
C-p,-»?) lCp,r;| = 1) and 

Af = -V'^ 2Gc Y [(^^i,'?J,rclp,^,fe,r') + (Cp,r,j,rC-p,,,,fc,r')) ^jkl^TT'Cp^r^ 

P,V,j,k,T,T' 

(6) 

To be precise, the index j in c^^j^ labels states of the 3 representation of sn(3), 

while in j ^ it labels states of the 3 representation. 

For convenience, we introduce the notation, m = (p,r/, r), fn = (— p, ry, — r). 
Clearly, fn = m. The BCS vacuum |$) is defined as the state which is annihilated 
by the operators djm such that 

dim = Clm - -fCm.(4m ~ 4m)' ^' < m, 

dirn = c\^ + Krn{c2m- c-m), m < 9.' , (7) 

that is, djm\^) = 0. The parameters Km, Km are real. The notation Q' < m 
means symbolically that m is such that \/p'^ + — > 0, while m < Q' means 
that ^/jPTIP - /i < 0. The expressions for d2m, d^m, are obtained by circular 
permutation of the indices 1,2,3. The transformation d?]) is not canonical, since 
{dim,dj^} 7^ Sij6mn, but the corresponding canonical transformation, which is not 
needed for the present purpose, may be easily obtained. 
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At this point, a short explanation may be in order, concerning the generahzed 
Bogohubov transformation, symmetric in the color indices, which has been proposed 
in eq. (25) of the first paper cited in [11], and which was supposed to "diagonal- 
ize" the pairing hamiltonian of the Bonn model. However, there is an error in that 
equation and the BCS vacuum associated with the Bogoliubov transformation it de- 
scribes fails to produce a non-vanishing gap A. That transformation is appropriate 
to diagonalize a pairing Hamiltonian which, although seemingly similar, is actually 
essentially different from the pairing Hamiltonian of the Bonn model, since that 
Hamiltonian is not invariant under SU{3). That is the reason for the corrigendum 
in [n]- 

The color symmetrical BCS state reads 

1$) = exp I Krr^A]^ + Krr^A.m] \0n'), (8) 

j=l \n'<m m<Q' J 



where 




|oo') = I IT IT 4m4^ 1 



and 

|0) denoting the quark vacuum. The expressions for ^2m'^3m' obtained by 
circular permutation of the indices 1,2,3. For simplicity, pairing operators involv- 
ing anti-quarks c (negative energy states) are not shown, but, in principle, their 
contribution should be included. In conformity, we will not show the contribution 
of anti-quarks to color superconductivity, but it is straightforward to include that 
contribution. 

We use the notation {W) = {ft\W\<^) / . We easily find 

(cUm) = -^-^^, z/,, (sWn) = 3^^, n'<m. (9) 
On the other hand, 

(4^.™) = -r^f^, i + J, {cUc,m) = 1 - T^^, rn < O'. (10) 
1 + 3Ki 1 + 3K^ 

The derivation of of eqs. ^ and (jlOp is summarized in the Appendix. 
The U (3) generators read 

Sij = ^c|m'^J"i- (11) 
m 

Clearly, the generators 5*^^ of SU{3) considered in the Introduction are related to the 
generators 5^;. For instance, Sx^ = S12 + S21, = -^(5'i2 -5*21), Sx-^ = 5'ii-S'22, 
etc. 

Since 

m<n' ^ m>n' ^ + -^^rn 

we may obviously insure that (Sij) = 0, i ^ j, which is the condition for color 
neutrality, by conveniently choosing K^, K^. 
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Next we compute the contractions {c2mCvm) = {c^mC2m) = (cimCsm) = {cTmCim) = 
{c3mC2m) = {cimCsm) =■ Dm, where Dm is real. We have, 

Km K 
1 + 1 + 3^2^ 

the derivation of these relations being left to the Appendix. 

We are now able to compute the expectation value of Km fa- We obtain 



+3 ^™ 1-^ 



2Kl \ ^ Km \ ,,3A 
+ 2A ^ ] + V- 



Ai = As = A3 



+ 3Kl J 1 + 3i^2^ J 4:Gc 

(14) 

^(ET^+Er^)-:A, (15) 



where stands for \/ p'^ + M'^ — fx. It is convenient to define the angles 6m, 6m such 



that sin0^ = ^Kml^/l + 3Kl, cosdm = l/yT+SKl, sm6m = V^Km/Jl + 3Kl, 



cos = 1/ Y 1 + 3ir^. Then, we have 

{Kmfa) = ^ (^Em sin^ 6*^ + 2V3A sin 6*^ cos dm^ 

(/ \ \ 3A2 

Em (3-2 sin^ ^^J + 2V3A sin 6m cos ^^j + V (16) 



VSA = I ^ sin 6*^ cos 0m, + ^ sin^mcos^mj • (17) 

\n'<m m<Q' J 

Having in mind eq. (fT2]) . the color neutrality constraint (Sij) = reduces to 

- sin^ 0m + ^ sin^ §m = 0. (18) 

Q'<m m<C 

The extremum condition reads (see Appendix) 

cos 26m = — ™ , sin20„, = ^ , 0' < m, 

V(em-A)2 + 3A2' v'(^m-A)2 + 3A2' 

o2 em - A ~ V3A , 

cos 20„, = , , sm 2dm = , , m < U , 

V(em-A)2 + 3A2' V(em-A)2 + 3A2' 

(19) 

where A is the Lagrange multiplier which ensures the color neutrality constraint 
(jlSp . The gap equation for color neutral superconductivity reads 



By setting A = 0, the gap equation for color neutral (not color symmetrical) 
superconductivity is obtained. Notice that we have been able to achieve color 
neutrality, eq. ([1]), independently of the A value, without introducing any extra 
Langrange multiplier, as opposed to what is done in O [10]. We stress, however, 
that the implementation of ()17p is essential. Otherwise, the present BCS vacuum 
can be reduced to the usual one, discussed in the next section, by an appropriate 
color rotation. 

At the extremum, the expectation value of Kmfa reduces to 
{Kmfa) = Yl (2^™ - ^m) + Yl -Em) + V 

m<Q' m>Q' ^ 

= 3Ye^m+Yl Yl (^"^ -Em)+V (21) 

m<n' m<n' n'<m ^ 



where Em = sj {£m — A)^ + SA^. In terms of proper canonical quasi particle opera- 
tors fim such that {firmfjn} = ^ij^mn, the constrained Hamiltonian reads 

Kmfa = ^ + ^(^m - Em) + V — ^ X] I Emfjmfjm + ^mflmfsm 

m<Q' rn ^ m \j=l 

(22) 

As explained in [11], the index j in fjm does not specify a well defined color. Indeed, 

flm = l^l{dlm — d2m), /2 m = l^2{dlm + C?2m — 2(i3m), /sm = ^3(^1^ + ^2™ + f^^Sm), 

where ki, K2; ^^Sj are normalization constants. 



3 Comparison with usual color superconduc- 
tivity 

In the usual approach to color superconductivity, which breaks color symmetry 
O [6], we have A3 = A 7^ 0, Ai = A2 = 0. The BCS transformation which 
digonalizes Kmfa reads 



with 



2 _ ^2 _ 2a 3 - ^ +6^ - 1 

The gap is expressed as 

A = 2^y2 amf3m, 

m 

SO that, the gap equation reads 

Gcy^^^ (23) 
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The constrained mean field Hamiltonian reduces to 



KmFA = 3 ^ em + X] ~ + ^{^m- Em) + V 

m<Q' m<Q' Q'Krn ^ 

m 

The gap A is only due to the pairing correlations between colors 1 and 2. The 
gap equation has essentially the same form in the color symmetric (eq. (j20p l and in 
the color symmetry breaking (eq. (j23p ) description of the color superconductivity. 
However, in the first case we have automatically (Su) = {S22) = {S33), while in the 
second case we have (^n) = {S22) / {Sss)- In [5, it has been argued that the usual 
color superconductivity, where only two colors participate in the gap, softens the 
equation of state. It is expected that the participation of the three colors on the 
same footing will have significant consequences for the high density phases of QCD 
which exhibit color superconductivity. We do not discuss the important color-flavor- 
locking mechanism, which also leads to color neutrality, since our main concern is 
the two flavor case. 




Figure 1: Lowest energy in the BCS approximation for the schematic Bonn model [TT] 
with Q = 10, versus the quark number, in arbitrary units. Thick line: the lowest energy 
of the color symmetric sector, according to the new method presented in Section 2; thin 
line: the result of the usual approach, which is summarized in Section 3, and describes 
the full groundstate energy of the model, including the "unphysical" sectors. 



4 Conclusions 

Quarks in QCD are, in the confined phase, only allowed to form colorless states. 
However, at high temperature and density, quarks are expected to be free. It has 
been shown that in the deconfined phase, which then prevails, the tendency for the 
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formation of BCS pairs occurs O |^. This gives rise to the so called color super- 
conducting phase, which, in the usual treatment jSj [H EJ (U [U [10], violates color 
symmetry. In this paper, we show that color superconductivity is not incompati- 
ble with color symmetry. It will be interesting to compare, in a realistic model of 
the NJL-type, the properties of the phase described by the color symmetric version 
of color superconductivity presently considered, with the corresponding properties 
described by conventional color superconductiviy, which drastically breaks color 
symmetry. This has been done for a schematic model in ref. [11], and, there, the 
difference found is important. In the Figure, we illustrate the performance of the 
new method here propose, and compare it with the usual approach to color super- 
conductivity, in the context of the QCD inspired schematic model considered in [llj . 
This model, which is characterized by a parameter Q measuring the level degener- 
acy, is admittedly unrealistic, but is quite useful to test approximation techniques, 
since it is analytically solvable. Applications to realistic situations are in progress. 

An effective QCD theory, as is appropriate to describe e.g. the interior of neutron 
stars, will have vanishing confining force at high temperatures and densities, due to 
asymptotic freedom, but should also be consistent with color singlet-ness [3 [8]. The 
phase of superconducting color symmetric states is supposed to exist in the interior 
of neutron stars with high density, where the simplified NJL model becomes identical 
to an effective QCD field theory and, thus, realistic. 

The present approach is in contrast with those followed in refs. [13] and [9], 
where the color neutrality problem was previously addressed. In [13], the authors 
resort to rather involved projection techniques to extract color neutral states out 
of BCS sates which violate color symmetry. It should be pointed out that the 
correlations described by the present approach need not coincide with those arising 
within the framework of the projection technique. In [9], color neutrality is defined 
by the condition that the average or expectation value of some of the eight Cell- 
Mann operators vanishes, that is, color neutrality is implemented with the help of 
appropriate Lagrange multipliers. In connection with this question, a reference to 
a recent work of Buballa and Shovkovy [10] is appropriate. These authors observe 
that, if a common chemical potential /i is used for all colors, then the quark numbers 
(S'li), {S22), (5*33), are not all the same, and hence will violate color neutrality. 
While, in the present approach, equality of the quark numbers for different colors 
is achieved with a common chemical potential, even if we set A = 0, in [10] one 
advocates using different chemical potentials for different colors, to insure (^n) = 
{S22) = {333)- However, in [10], it is already observed that this condition is not 
sufficient to insure color singlet-ness, due to its instability under color rotations. We 
have proposed a solution to the emerging problem. 

The present article also provides theoretical "tools" for constructing phases with 
superconducting color-symmetric states. 

Although our aim is not to discuss the QCD Meissner effect, we observe that 
this effect may be treated in terms of a Hamiltonian of the form 

H = jd^x Q(B, • + E, • E,) + n* • n + (DA)* • (DA) + V{A* ■ A)^ , 

describing the color superconducting phase in interaction with the gluon field. In 
the definition of the covariant derivative D one should keep in mind that A belongs 
to the 3 representation. It is clear that the standard treatment goes through with 
our approach. We observe that our color symmetric BCS theory should lead to 
qualitatively the same masses for the gauge bosons as are found in the literature. 
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The only place where an effective model such as the NJL model comes in is in the 
estimation of the function V which depends on | Aip, IA2P, IA3P. Our approach 
ensures that V will depend only on the combination |Aip + IA2P + IA3P. 



Appendix 



Derivation of eqs. ([9]) and (1101). 

For W < m, we have 




implying 



= -Kl + {N^ -X^), N^ = 2Kl - 2KI (iV„ - , 

which leads to X„, = - K^J {1 + 3K^,), = '2,K'^/{1 + 3if^). The corresponding 
expressions for m < Q' are similarly obtained. 



Derivation of eq. (1131). 

For W < m we have 



3m/ / ' 



which imply 



-Dm — Km '^KmDm ~\~ KmP-mi Pm — "^KmPmi 



where Pm = (cimCim) = (c2mC2rn) = (csmCsm), is also real. The procedure for 
$7' < m, is analogous. Finally, we find Pm = and eq. ()13p follows. 



Proof of eq. (fT9l). 

Let 

^{e, e, A, A) = {Kmfa) + 2A - ^ sin^ 9m. + ^ sin^ 9m. , 

\ Q'<m m<n' J 

where {Kmfa) stands for its expression in eq. p6|) and A is a Lagrange multiplier. 
We find 

Se,^^ = 4(em - A) sin 9m cos 9m + 2\/3A(cos^ 9m - sin^ 9m) = 0, Q' < m, 
dgjl! = -4:{£m - A) sin 9m cos 9m + 2\/3A(cos^ 9m. - sin^ 9m) = 0, m < Q'. 

The condition (9a^ = yields eq. (fT7|) . and the desired result, eq. (fTUj) . follows now 
easily. 
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